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Abstract 

A new physical implementation for quantum computation is proposed. The vibrational modes 
of molecules are used to encode qubit systems. Global quantum logic gates are realized using 
shaped femtosecond laser pulses which are calculated applying optimal control theory. The scaling 
of the system is favorable, sources for decoherence can be eliminated. A complete set of one and 
two quantum gates is presented for a specific molecule. Detailed analysis regarding experimental 
realization shows that the structural resolution of today's pulse shapers is easily sufficient for pulse 
formation. 
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The possibilites offered through quantum computation have been well known for almost 
twenty years [1, 2]. In the last few years the first basic ideas for realization have been 
discussed and experimental efforts have been made using cavity quantum electrodynamics 
[3], trapped ions [4, 5], and nuclear magnetic resonance [6, 7]. 

In this Letter we demonstrate how to use vibrational modes of molecules in an ensemble 
for complete quantum computation processes. In an N-atomic molecule 3N-5 (or 3N-6) 
normal modes can be identified, and each of these normal modes can be used to define 
a quantum bit. Two different excitations in each mode are referred to as |0) and |1). To 
implement quantum logic operations in this qubit system, shaped femtosecond laser pulses in 
the IR-regime are used. They are designed with the help of optimal control theory (OCT) 
[8, 9] which are in principle realizable in pulse shaping experiments [10, 11]. Interaction 
between different qubits is an inherent property of the whole qubit system belonging to 
the same molecule. For designing globally applicable laser pulses in just one qubit, this 
interaction has to be suppressed; in the case of conditional gates this influence is explicitely 
used. Both cases can be realized in our OCT-functional. 

For initial state preparation again shaped femtosecond laser pulses are used. They can 
be calculated with our standard OCT-algorithm which can also compensate for rotations of 
molecules in the gas phase. [12, 13]. Vibrational state preparation is optionally possible for 
the vibrationless ground state, fundamentals, overtones and combination modes. In our work 
the initial state preparation is a preparation of a combination mode for which maximally 
two shaped femtosecond laser pulses of different wavelengths are needed. 

The detection of states after a computation can be achieved by applying standard laser 
diagnostics (for example time resolved IR-spectroscopy [14], laser induced fluorescence spec- 
troscopy [LIF] [15], stimulated emission pumping [SEP] or resonant multi photon ionization 
[REMPI] [16]). 

In the following part of the Letter we will show that vibrationally excited molecules 
provide a realistic physical system to implement a quantum computer. We will demonstrate 
a realization of the above ideas using acetylene {C2H2) as a two-qubit model system. Our 
new OCT-functional, laser pulses for elementary global one and two quantum gates (which 
form a universal quantum gate altogether [17]), and the preparation of entanglement will be 
presented. An analysis of the shaped femtosecond laser pulses will show that these pulses 
are simply structured. The problem of decoherence is discussed and a solution is found. 
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Error detection is also feasible. 

OCT is a very powerful tool for calculating laser pulses guiding a quantum system to 
any selected objective, in particular the preparation of eigenstates is enabled. In this case 
the goal in OCT is to find a laser field e(t), which drives a system from initial states 
V'ifc(O) = 4>ik a t time t — to final target states <pfk at a fixed time t = T. The initial 
and target states correspond with the eigenstates defining our qubits. The calculated laser 
pulses e(t) represent the quantum gates. 

Algorithmic schemes allow the formulation of the optimization problem in terms of the 
maximization of the functional 



The first term gives the overlap between the laser driven wavefunction ipik(t) and the 
desired target states <f>f k ; the second one represents the laser field e(t) which drives the 
system wavefunctions towards the target states. The laser field intensity is limited by the 
penalty factor a. The last term ensures that the time dependent Schrodinger equation is 
fulfilled. H Q and the dipole moment vector field ji have been calculated ab initio, which 
means that all molecular characteristics of the electronic and nuclear movement and their 
interactions are exactly implemented in the calculations. 

Varying the above functional with respect to the variables ipik{t), ipfkif), and e(t) and 
searching for 8K = leads to a set of 2k+l coupled differential equations. The solution is 
found iteratively and results in the optimal quantum gate e(t) as a self consistent solution 
to this system. 

The new feature which is of highest importance in the context of quantum computation 
is the fact that we are able to optimize k transitions simultaneously with the above new 
functional. Its recent implementation in our algorithms enables us to calculate globally 
applicable quantum gates, because for each basis state of our qubit system the correct tran- 
sition k is optimized, and undesired transitions can be explicitely suppressed. For example 
in a two-qubit system a global NOT-pulse in the second qubit has to fulfill 1 00) — > |01), 
|01) -> |00), 1 10) -> |11) and |11) -> |10) at the same time. 

These basic principles for concrete quantum computation with vibrational modes are now 
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discussed for an example system: acetylene. C2H2 contains five normal modes (ni^n^n^^) 
in the electronic ground state (n 4 and n 5 are doubly degenerate). The asymmetric CH- 
stretching (n 3 ) and the ds-bending mode (n 5 ) are IR-active. We have chosen them to define 
our two-qubit system [18] because direct experimental access is possible. One quantum in 
each mode is referred to as |0), two quanta as |1) (other classifications would also have been 
possible). Our complete two-qubit system is described with the following four basis states: 
|00) := (00101); |10) := (00201); |01) := (00102); |11) := (00202). 

We have computed a two dimensional ab initio potential energy surface (PES) using a 
standard quantum chemical package (GAUSSIAN98 [19]). One system-adapted coordinate is 
sufficient for the description of one IR-active mode, because the relative distance between the 
two hydrogen atoms and the relative distance between the two carbon atoms rests constant 
during the motion (coordinates: R — > ds-bending mode, d — > asymmetric CH-stretching 
mode, see [12]). In order to calculate the vibrational eigenstates (j>f k we solved the nuclear 
Schrodinger equation for our PES on a grid applying a relaxation method followed by an 
additional diagonalization step [20, 21]. We calculated 80 eigenfunctions up to an energy of 
23150 cm -1 (the minimum energy is set to zero), four of them define our qubit basis. For the 
laser-matter-interaction we have computed the ab initio dipole moment vector field, which 
can be decomposed into two parts, each part being responsible for the excitation of one 
IR-active mode. This decomposition and the energetic spacing between the two vibrations 
ensure a very specific excitation of each mode (ds-bending mode 727 cm -1 (13.75 /xm), 
asymmetric CH-stretching mode 3289 cm -1 (3.04 um)). For the construction of our quantum 
gate e(t) the set of coupled differential equations is solved quantum dynamically applying 
the split operator scheme [22] for the wave packet propagations. 

With our new functional we have constructed global laser pulses for the following ele- 
mentary one-quantum gates: identity operation 1 (trivial case or explicitely computed), bit 
flips (NOT), phase shifts of tt and the Hadamard transformation in the second qubit (the 
case for the first qubit is analogous [23]): 

1: |0)(0| + |1)(1| 
NOT: 1 0> < 1 1 + |1)(0| 
IT: |0)(0|-|1)(1| 
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Hadamard H\0) -> -^= [|0) + |1)] and -> -L [|0) - |1)] 

v 2 v 2 

As a conditional two-quantum gate we have calculated a controlled NOT (CNOT) with the 

first qubit acting as the control state: 
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Exemplarily the optimized laser field for the CNOT gate is presented in Fig. 1. The 
first part of the figure displays the quantum gate e(t) versus time, the second one depicts a 
fast Fourier transform to analyze the spectral width and components of the laser pulse. It 
is not Fourier limited and pulse shaping is needed to form this electric field with a broad 
spectral range. The last part of Fig. 1 analyzes the pulse structure in the time and frequency 
domain simultaneously mirroring the temporal ordering of the frequency components. The 
substructure is significant and rather symmetric in time. 

For the realization of all our quantum gates our method needs only one single shaped laser 
pulse; a preceding special state preparation or pulse sequence is not necessary as in other 
approaches like cavity QED, trapped ions and also NMR experiments. We have analyzed the 
mechanism according to which our quantum gates work. In all cases not only states belonging 
to our qubit system, but temporarily energetically higher and lower lying vibrational levels 
are populated following a ladder climbing scheme. An intermediate coherent state is built 
up before all population is transferred into the desired eigenstates. These eigenstates are 
stable, no refocusing pulses are needed as in NMR quantum computing. 

The ratios r (defined as \ (ip ik {T)\4> fk )\ 2 ) we can reach are: r(NOT) > 91%; r(n) > 97%; 
r(Hadamard) > 90%; r(CNOT) > 90%. (For comparison, in cases of state prepa- 
ration with k = 1 we can always reach efficiencies of > 98%.) In each case the 
laser pulse found by the OCT-algorithm consists of several structural units, and the 
total laser matter interaction time is about 700 fs. The range of the frequency 
width is between 200 cur 1 (7r-pulse) and 560 cm" 1 (CNOT-pulse). The strength 
of the electric field e is moderate except for the NOT-operation, but still below the 
ionization limit (maximum strength: e(NOT) = 0.0224 a.u. (7.146-10 13 W/cm 2 ); 
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e ( 7r ) = 0.0057 a.u. (1.1597-10 12 W/cm 2 ); e(Hadamard) = 0.0095 a.u.(3.226-10 12 W/cm 2 ); 
e(CNOT) = 0.0070 a.u.(1.749-10 12 W/cm 2 )). First investigations in other systems show 
that the field strengths can be lowered down [24, 25]. 

With our OCT-algorithm we are also able to calculate pulses which create entangled 
states. The efficiency or fidelity of our Bell state preparation (starting at one of our basis 
states and proceeding via a Hadamard-state) is r > 96% in each case. One part of the 
preparation of -j= [|00) + |11)] is shown in Fig. 2. 

The prepared Bell state is not a long-lived normal mode, but a superposition state. In that 
context it is relevant to investigate the time evolution of the corresponding wave function. 
In Fig. 2 and 3 we have depicted the time evolution of the real part of our wave packets after 
the preparation of a Bell state and a Hadamard state (of course, we have experimental access 
only to |"0| 2 ). The wave functions (real and imaginary part) develop in time systematically, 
and \ijj\ 2 comes back in predictable intervals. Their evolution in space is confined to a certain 
area. Once again one does not need refocusing pulses. 

To realize our calculated gates, pulse shaping techniques are needed. In the spectral 
range from 430 nm to 1.6 fim liquid crystal modulators (LCM) are already commercially 
available, having typically 128 discrete pixels. Powerful lasers operating in the mid-infrared 
already exist (free-electron laser FELIX (4-250 /j,m) [26]). Pulse shapers that work in the 
mid-IR and far-IR spectral ranges are actually under development; the development of new 
types of nonlinear crystals already shows promising results in those directions [27]. Anyhow, 
a realization of a pulse shaper in the mid-IR and far-IR spectral ranges will follow a pixel 
scheme. The mask function applied to the pulse shaper in order to generate the shaped laser 
pulses is the direct interface between theory and experiment. Therefore we calculated the 
mask functions needed to realize our quantum gates and investigated their complexity. The 
structure of the mask function (transmission and phase) is very simple and only a few pixels 
are needed. As an example the mask function for the CNOT is displayed in Fig. 4. The 
original pulse to be sent through the pulse shaper (Fig. 4 top) is a 171 fs strictly Fourier 
limited pulse (center frequency at 674 cm -1 , FWHM 590 cm -1 , maximum electric field 
strength 0.0280 a.u. (2.799-10 13 W/cm 2 )). 

The ratios we can achieve with our gates at this stage are rather high, but still below 
100%. We have analyzed the losses. Most result from leakage into vibrations that lie 
outside our chosen Hilbert space and therefore bear little disturbing influence for further 
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calculations. The effect is a diminishing signal when detecting results after a quantum 
computation process. An estimation shows that at room temperature, normal pressure and 
an assumed laser focus zone of 200 \im in diameter and 1 cm in length a few hundred 
quantum gates could be applied successively. This is already enough for some very useful 
devices (e.g. a quantum factoring engine). We have found that the losses are localized 
in very few specific vibrational states and can be detected there without disturbing the 
information stored in our qubit system. Thus error detection is possible in our proposed 
quantum computing model. 

The context of leakage is closely connected to the topic of decoherence. One could think 
of two different sources of decoherence: resonances between different vibrations (anharmonic 
resonances and Coriolis coupling) and collisions between molecules. Recent investigations 
have shown that anharmonic resonances do not cause decoherence at all, but sometimes 
provide easier mechanisms for designing quantum gates [24, 25]. Regarding molecules in the 
gas phase, the number of collisions can be kept low, and typical lifetimes of vibrations are 
in the nanosecond regime or longer, depending on the experimental conditions. The more 
significant problem is motion of the molecules within the laser-matter-interaction volume 
and possibly out of that zone before the next laser pulse is applied. The repetition rate of the 
lasers ist too small and diminishes the number of possible steps in a quantum computation 
process significantly. These problems can be avoided if fixed molecules are used. One can 
also think of specially designed macro-molecules adsorbed on a surface with the active centers 
still free to vibrate [28, 29, 30] or of molecules embedded in films [31]. Molecular design 
techniques for self assembling molecules are well developed at the frontier of nano technology 
[32, 33]. Rotations are also suppressed in fixed molecules, and therefore decoherence due to 
Coriolis coupling is precluded. 

In the future, the design of macro-molecules consisting of repeated subunits, e.g. sub- 
systems connected through conjugated CC-bonds, might play an important role for molec- 
ular quantum computing processes. The interaction within each subunit and between the 
subunits is secured using the system-inherent influence of the vibrational modes upon one 
another as explained above. In particular, a collective vibrational mode of all subunits could 
be used as a data bus in such a quantum network. Every subunit and each vibrational mode 
has access to the same information coded in the excitation of this collective vibration. An 
N-atomic macro-molecule is an immense qubit system: the number of utilizable qubits scales 
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with 3N, in NMR and ion trap experiments with N. 

In summary, vibrationally excited molecules are well suited to implement quantum com- 
putation processes. Shaped femtosecond laser pulses act as extremely fast quantum gates 
and can be designed by applying optimal conrol theory. The calculated laser pulses are 
simply structured and could be realized in an experimental surrounding where decoherence 
is negligible. 
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Figure captions 



FIG. 1: Controlled NOT gate for a two-qubit system in vibrationally excited acetylene. Depicted 
are the strength of the electric field, a fast Fourier transformation, and a FROG-representation. 
The efficiency of the gate is > 90%. 



FIG. 2: Preparation of a Bell state: ^ [|00> + |10)] [|00> + |11)], efficiency > 98%. R is the 

coordinate of the cis-bending mode, d describes the asymmetric CH-stretching mode. Displayed 
are the laser pulse and the time evolution of the real part of the wave function during the pulse 
and afterwards. The * indicates the initial and target state, the box labels states with the same 



FIG. 3: Preparation of a Hadamard state: |00) -> [|00) + |01>] , efficiency > 99%. R is the 
coordinate of the cis-bending mode, d describes the asymmetric CH-stretching mode. Displayed 
are the laser pulse and the time evolution of the real part of the wave function during the pulse 
and afterwards. The * indicates the initial and target state, the box labels states with the same 



FIG. 4: Mask function for the CNOT-gate. For explanation see text. 
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